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Abstract. 

Josephson weak links in superconductors can be engineered such that the phase 
difference across the junction is it (jr-j unction), leading to the reversal of the cur- 
rent. The conditions under which the analogous effect of supercurrent counterflow can 
be achieved in a double-well atomic Bose-Einstein condensate are investigated. It is 
shown that this effect is observable for condensates of up to a few thousand atoms, 
which are initially prepared in an anti-symmetric '7r-state' and subsequently subjected 
to a uniformly increasing magnetic field gradient. This effect is found to be only 
weakly-dependent on trap geometry, and can be observed in both attractive and re- 
pulsive condensates. 



Note: A substantially revised version of this manuscript (with various 
figures removed, improved discussion and some additional material) is 
available on cond-mat/0312396 
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1. Introduction 

The creation of superconducting P and superfluid [2j weak links has led to the 
experimental observation of Josephson effects [3], arising as a result of macroscopic 
quantum phase coherence. Josephson weak links are typically created by connecting 
two initially independent systems (superconductors / superfluids) via a barrier with 
dimensions of the order of the system healing length. Such junctions lead to a 
variety of interesting phenomena [I], such as dc- and ac- Josephson effects, Shapiro 
resonances, macrosopic quantum self-trapping and 7r-phase oscillations. Observations 
in superconductors preceded those in superfluids, due to the much larger healing 
lengths, thus enabling easier fabrication of weak links. Evidence for Josephson-like 
effects has been observed in 4 He weak links f5J, and unequivocally demonstrated for 
weakly-coupled 3 He systems jO]. The recent achievement of dilute trapped atomic Bose- 
Einstein condensation (BEC) [7j gives rise to a new system for studying Josephson 
effects. In particular, such systems enable the investigation of dynamical regimes 
not easily accessible with other superconducting or superfluid systems. The simplest 
atomic Josephson junction can be realized by a condensate confined in a double-well 
potential. To allow control of the tunnelling rate, such a system can be constructed by 
raising a barrier within a harmonic trap containing an atomic condensate; this can be 
achieved by the application of a far-off-resonant blue-detuned laser beam, which induces 
a repulsive gaussian barrier jHj. Atomic interferometry based on such a set-up was 
recently reported [Oj. Alternatively, a condensate can be created directly in a magnetic 
double- well structure JU]- Remarkable experimental progress has led to the creation 
of atomic BEC Josephson junction arrays, in which the harmonically trapped atoms 
are additionally confined by an optical lattice potential, generated by far-detuned laser 
beams. Phase coherence in different wells was observed by interference experiments 
of condensates released from the lattice [TT| . In addition, Josephson effects [12] and 
the control of tunnelling rate has been demonstrated ^J^l- Although experiments 
(and theoretical analysis) of such systems are well underway, deeper insight into the 
diverse range of Josephson phenomena can be obtained by looking at the simplest single 
junction, double-well system. This system has already received considerable theoretical 
attention, with treatments based on a two-state approximation [T5 t lTo ^ lTT t lT ^ IT ^ l2T) t l2T] . 
zero temperature mean field theory [22123123122112111211 , quantum phase models [2~Sll2Ti] 
and instanton methods [HO] . 

This paper focuses on a double-well atomic BEC, and investigates the conditions 
under which the Josephson current can be engineered to flow in a direction opposite 
to that minimizing the potential energy. This phenomenon bears close analogies to 
superconducting 7r-j unctions [3X], in which a macroscopic phase difference <p = n is 
maintained across the superconducting weak link. Such behaviour has been observed in 
a variety of systems, as a consequence of different microscopic mechanisms. For example, 
7r-junctions in ceramic superconductors have their origin in the symmetry of the order 
parameter [32], with their experimental detection being central to the understanding 
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of high-T c superconductivity, it -junctions can also be created in ferromagnetic weak 
links [33], or by magnetic impurities |2]. Recent interest has focused on the creation 
of controllable 7r-junctions in superconducting/normal- metal/superconducting links, in 
which the current direction can be reversed by the application of suitably large voltage 
across the link 12^3- The reversal originates from the fact that the addition of an extra 
phase factor 7r is equivalent to reversing the sign of the current J c , since the superflow 
current obeys / = / c sin(0). A controllable 7T-SQUID (Superconducting Quantum 
Interference Device) has been demonstated [HBj, and it is clear that the manipulation 
of multiple such 7r-junctions will be important in the domain of quantum electronics. 
For example, an array of alternating — ir junctions allows the spontaneous generation 
of half- integer flux quanta. Such a circuit of multiple successive 0-tt junctions has 
been recently created between thin films of high-T c and low-T c superconductors |37j . 
generating a one- dimensional array of Josephson contacts with alternating signs of 
current. The superfluid analogue of a supeconducting 7r-junction is a metastasble ir- 
state. This was recently observed in 3 He weak links |HH] upon exciting the system by an 
oscillating driving force. Atomic BEC junctions behave similarly to those of 3 He. For 
example, in the usual manner of considering mechanical analogs of Josephson junctions, 
superconducting Josephson junctions can be mapped onto a rigid pendulum, whereas 
atomic tunnel junctions ( 3 He, BEC) map onto a non-rigid pendulum [SH], thus exhibiting 
richer oscillation modes. This model has been used to discuss so-called n oscillations 40 j 
and their stability in atomic BEC's while such states have also been shown to arise 
within the framework of an exact quantum phase model 

In this paper, we investigate the circumstances under which one can reverse the 
direction of the atomic current across a suitably-prepared condensate-condensate weak 
link, by the application of a linear potential gradient. We find that a BEC confined in a 
double well configuration can, for small values of the potential gradient, move towards 
the higher potential well, a phenomenon henceforth referred to as Josephson counterflow. 
Although such Josephson counterflow bears close analogies to the behaviour observed 
in superconducting weak links, we should point out that the counterflow discussed in 
this paper is 'global' (i.e. flow of entire quantum gas in opposite direction), as opposed 
to 'local' counterflow (i.e. across a single junction) in a superconducting 0-n junction. 

This paper discusses in detail the phenomenon of Josephson counterflow for the 
lowest state exhibiting such behaviour, namely the anti-symmetric first excited 7r- state, 
which is most amenable to experimental observation. Atomic counterflow dynamics 
are investigated in terms of experimentally relevant parameters, such as interaction 
strength, harmonic trap aspect ratio and gaussian barrier geometry. Our analysis is 
based on numerical simulations of the Gross-Pitaevskii equation in three dimensions 
and leads to the conclusion that there exists a realistic window of parameters in which 
atomic Josephson counterflow can be experimentally observed. This effect is found to 
be only weakly dependent on the strength of the interactions and can, in prinicple, be 
observed for both attractive and repulsive condensates. One should note the distinction 
between the 7r-state considered in this paper which is a time-independent solution, as 
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opposed to the it oscillations which arise as a result of a superposition of ground and 
first excited states. In an experimental realization, it may be difficult to create a pure 
7r-state, and the system may experience a combination of counterflow and 7r-oscillations. 
In this paper we show that, by careful initial state preparation, one can decouple the 
timescales for these two effects, and even induce counterflow in a direction opposite to 
that of 7r-oscillations, thus demonstrating the different origin of these two phenomena. 
Note that Josephson counterflow has already been predicted in condensates trapped in 
optical lattices, as a result of the renormalization of the mass in the lattice, based on 
Bloch wave analysis |42j . Such an interpretation is, however, not easily transferable to 
the double-well system. 

This paper is structured as follows. Sec. 2 introduces our main formalism, 
outlining the behaviour of a double- well condensate. Sec. 3 discusses the dynamics 
of atomic Josephson junctions, whereas the possibility of experimental observation of 
this phenomenon in current BEC set-ups is analysed in Sec. 4; finally, we conclude in 
Sec. 5. 



2. BEC in a double-well potential 

At low temperatures, the behaviour of a Bose-Einstein condensate is accurately 
described by a nonlinear Schrodinger equation known as the Gross-Pitaevskii (GP) 
equation. Throughout this paper we work in dimensionless (harmonic oscillator) 
units, by applying the following scalings: space coordinates transform according to 
r\ = a±~ 1 r i , time t' = u±t, condensate wavefunction ip'(r',t') = ^J~aJ^ip(r,t) and 

energy E' = ( y hu±)~~ 1 E. Here a± = ^Jh/mu>± is the harmonic oscillator length in the 
transverse direction(s), where u>± the corresponding trapping frequency. We thus obtain 
the following dimensionless GP equation (primes henceforth neglected for convenience) 
describing the evolution of the condensate wavefunction (normalised to unity) 

-l:V 2 + V(r) + ~g\ij(r,t)\ 2] 



id t i(j (r, t) 



i>(r,t). (1) 



2 

The atom-atom interaction is parametrized by g = g/(a\luo±), where g = J\f(4:7ih 2 a/m) 
is the usual three-dimensional scattering amplitude, defined in terms of the s-wave 
scattering length a, and M is the total number of atoms (mass m) . V (r) represents the 
total confining potential. Steady state solutions of the GP equations can be obtained 
by substituting ip (r, t) = e' t>Mt (f) (r). 

In the double-well configuration, the total confining potential is given by (see Fig. 
1(a)) 

V(r) = \ [(x 2 + y 2 ) + \ 2 z 2 } + h exp [-(z/w) 2 } + Az . (2) 

The first term describes a cylindrically symmetric harmonic trapping potential, with a 
trap aspect ratio A = uj z /lj±: the trap is spherical for A = 1, 'cigar-shaped' for A < 1 
and 'pancake-like' for A > 1. The second term describes a gaussian potential of height 
h generated by a blue detuned light sheet of beam waist w in the z direction, located at 




Figure 1. Double well potential with corresponding eigenenergies and eigenstates. 
(a) Schematic geometry of the total confining potential in the axial direction for a 
Gaussian barrier (height h = 4huj±, waist w = a±) located at the centre of the trap. 
Plotted are the symmetric (A = 0, solid line) and asymmetric (A = 0.5(%u)±/a±), 
dashed line) case, (b)-(c) Corresponding eigenenergies and eigenstates for the double- 
well as a function of the potential gradient A: (i) ground state (lower thin solid line), 
(ii) anti-symmetric first-excited 7r-state (thick solid line), (iii) first excited state with 
unequal populations, having more population in left well (dotted), or in right well 
(dashed). Parameters used here are g = ir and spherical trap geometry (A = 1). 



z = 0. Finally, the contribution Az corresponds to a linear magnetic field gradient A 
pivoted at the centre of the trap. The populations of the two wells are equal for A = 
(symmetric case, solid line in Fig. 1(a)), whereas A > (dashed line in Fig. 1(a)) 
leads to a tilted potential, which induces tunnelling. The application of the magentic 
field gradient A > additionally shifts the trap centre to the z > region. However, 
this shift is negligible for the parameters studied throughout this work, and will be 
henceforth ignored. 

The eigenenergy curves of the double-well condensate are calculated by numerical 
solution of the time-independent GP equation, as discussed in more detail in our 
preceeding paper j2Zj. Sufficiently large interactions lead to the appearance of a loop 
structure The loop structure for the first excited state is shown in Fig. 1(b). 

Corresponding wavefunctions for ground (solid) and first excited state (dashed) are 
shown in Fig. 1(c) for A = 0. The three eigenstates are (i) a symmetric state with 
equal population in both wells (solid line), (ii) an anti-symmetric V-state' with equal 
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population in both wells and a phase difference of it across the trap centre, and (iii) 
two higher energy 'self-trapped' states with most of the population in either the left 
(dashed) or the right (dotted) well. 

In this paper, we are mostly interested in the behaviour of the 7r-state. We will 
show that, if the system is prepared in the n state, the subsequent temporal evolution 
of the system can exhibit Josephson counterflow. 

3. Josephson Counterflow Dynamics 

To demonstrate counterflow, we study Josephson dynamics of the 7r-state under a 
potential gradient A > 0, which is linearly increasing at rate R (i.e. A = Rt), such 
that the right well lies higher than the left well (dashed line in Fig. 1(a)). On the 
basis of the usual Josephson relations, one might naively expect the atomic current 
to flow towards the region of lower potential, i.e. towards the left well. Instead, we 
observe that, for suitable parameters (see later), superflow can occur from the lower 
potential energy well (left) to the higher potential energy one (right). This is a direct 
consequence of the n phase difference of the antisymmetric initial wavefunction, and 
does not occur for a system in the symmetric ground state (for which flow always occurs 




A r r 

Figure 2. (a) Evolution of fractional population difference N as a function of potential 
gradient A, for a system initially prepared in a 7r-state. Here h = 4?ilo±, X = 1, g = it 
and the potential gradient A = Rt increases at constant rate R = 10~ 3 (hui'^/a±). The 
corresponding population differences for the ground state (thin solid), first excited 
(dashed) and second excited state (dotted) eigenstates are also shown, (b) Snapshots 
of the evolution of the density distribution for case (a) when (i) A = 0, (ii) 0.1, (iii) 
0.3 and (iv) 0.8 in units of (hiv±/a±). The population of both wells is initially equal 
(A = 0). As the gradient is increased in (ii), (iii), population starts being transferred 
towards the right (z > 0), upper well. Increasing the asymmetry beyond a threshold 
value leads the population to be once again transferred to the left (z < 0), lower well. 
Eventually, (iv), a transition to the second excited state occurs (d) (see movie). 



Josephson Current Counterflow in Atomic BEC's 7 




-0.3 0.3 -0.3 0.3 -0.3 0.3 

AAA 

Figure 3. Evolution of fractional population difference N as a function of potential 
gradient A for identical trap configurations (A = 1, h = Ahuj± and R — 10~ 3 (hui^/ a ±)) 
and increasing nonlinearity (a) g = tt, (b) 4tt and (c) 107r. Solid bold (dashed) lines 
indicate corresponding eigenstate populations for the ground (first excited) state. 

towards regions of lower potential energy - see, e.g. |2Zj)- Our study focuses on the 
dynamics of the fractional relative population, N = (Nl — Nr) / (Nl + Nr), between the 
two wells, as opposed to the current through the junction (which is the derivative of N). 
The dependence of N on potential gradient A is shown for a spherical trap in Fig. 2(a), 
for a system initially prepared in the 7r-state. In such a state, the population of both 
wells is initially equal (A = 0). As the gradient is increased, population starts being 
transferred towards the right, upper well, such that we observe Josephson counterflow to 
regions of higher potential energy. However, increasing the gradient beyond a threshold 
value leads to suppression of this effect, with the population once again transferred to 
the left, lower well. Eventually, the perturbation due to the applied potential becomes 
so pronounced, that a transition to the second excited state occurs [2Z1- Characteristic 
density snapshots of this evolution are shown in Fig. 2(b). The initial counterflow 
dynamics can be understood by means of lowest order perturbation theory. However, 
such a simple picture no longer gives an accurate description when the population 
difference becomes large. The two-state model [1311111121111111111201 also reproduces 
the initial counterflow dynamics. However, since this model contains no mechanism 
for removing the system from the counterflow state, the two-state model predicts 
equilibration in a macroscopically quantum trapped state with larger population in 
the higher well. This inadequacy of the two-state model is based on the fact that it does 
not take higher lying modes into consideration [2*7j . 

Increasing the nonlinearity causes a reduction in the amount of initial counterflow, 
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and thus tends to inhibit the experimental observation. Fig. 3 illustrates the reduction 
in counterflow due to increased nonlinearity for fixed trap geometry, with g increasing 
by a factor of 10 from (a) to (c). This would appear to restrict the observation of the 
phenomenon to moderate nonlinear it ies (see next section for experimental estimates). 
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Figure 4. Dependence of fractional population difference dynamics on trap 
geometry. Plotted is the evolution of N as a function of potential gradient A for 
R = 10~ 3 (tiu)j_/ a j_) (black), g = 107r and different trap aspect ratios: (a) cylindrically- 
symmetric trap A = 1 (h = 4Tlui±) (same as Fig. 3(c)), (b)-(c) 'pancake'-like traps with 
A = \/2 (h = 6Huj±) and (c) A = y/B (h = 15Huj±), respectively. Note that larger barrier 
heights have been used with increasing aspect ratio, such that the density minimum 
at the trap centre is roughly constant from (a)-(c). 

It is thus natural to ask if other factors (e.g. modifying initial trap aspect ratio, or 
changing barrier height h) will have the opposite effect on the amount of counterflow, 
hence enabling observation of Josephson counterflow even for large nonlinearities. For 
example, tunnelling has been predicted to be enhanced for 'pancake' traps (A > 1) [24J. 
Indeed, for weak nonlinearities, such traps lead to a slightly increased counterflow 
amplitude. Furthermore, such geometries feature enhanced energy splitting between 
ground and first excited state, thus making them more robust to coupling due to external 
(e.g. thermal [THlUnillEllH]) perturbations. However, the reduction of the amplitude 
of counterflow due to increased nonlinearities tends to largely suppress this geometry 
dependence, as shown in Fig. 4 for trap aspect ratios in the range 1 to a/8. In plotting 
this figure, the barrier height h has been increased for larger A, such that the peak 
density in the trap centre remains essentially unchanged. 

We should further comment on the extent to which our above findings depend on 
the rate R with which the linear potential gradient A = Rt is applied. The thick solid 




Figure 5. Evolution of fractional population difference N as a function of time 
(thick black lines) for different rates R of increase of the potential gradient: (a) 
R = 1CT 3 {hw\f a j_) and (b) R = 10~ 2 (hui^_/a±). Grey lines plot corresponding 
evolution of the population difference for the case when the potential gradient is held 
constant after a time t — (a) fOOwJ 1 and (b) lOwJ 1 Other parameters used, as in Fig. 
3(a). 



lines in Fig. 5 correspond to the evolution of the fractional population difference as a 
function of time. Fig. 5(a) shows the dependence for R = lO~ 3 (huj? L /a±) (as used in all 
earlier figures), whereas Fig. 5(b) shows the corresponding behaviour when the gradient 
is increased at a rate ten times faster than (a). One observes the following effects: 
firstly, the amount of maximum (initial) counterflow is significantly reduced (roughly 
by a factor of 2) by increasing the rate R by a factor of 10. Secondly, counterflow can 
only be observed for much shorter times (roughly reduced by a factor of 10). 

A final question of interest is what would happen to the population difference if 
the applied linear potential is ramped up to a particular value and subsequently kept 
constant. The most striking behaviour will occur when the gradient is kept constant at 
the point of maximum population difference, as indicated by the grey lines in Fig. 5. 
We see that in this case, the population remains trapped in the right upper well, i.e. 
macroscopic quantum self-trapping occurs to a state with higher potential energy. In 
this regime, where the gradient does not exceed the value at which the flow is reversed, 
the two-state model predicts the behaviour correctly. 

4. Experimental Considerations 

Having demonstrated the existence of Josephson counterflow for a 7r-state initial 
wavefunction, we now discuss the feasibility of such observation in atomic BEC 
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experiments. Firstly, we need to discuss how such states with a node in the wavefunction 
and odd parity behaviour can be generated. Although not necessarily the most efficient 
method, here we consider phase imprinting |44j . Starting from the condensate ground 
state in a harmonic trap, population can be transferred to the first excited state by 
applying a light-induced potential of the form 

V r (z, t) = a sin (7rt/r ) tanh (z) (3) 

for (t < tq) , where a and tq are constants which we vary. This is equivalent to applying 
a 7i phase shift to the system. At t = r , the potential V r is suddenly switched off, and 
the system exhibits Rabi oscillations of variable amplitudes and frequencies between the 
initial (ground) state and the first excited state. Fig. 6 shows fractional occupations 
((a), (c)) and corresponding population differences ((b), (d)) for two different initial 
state preparation cycles. The fractional occupation of state i, denoting here ground 
(upper, thick solid) or first excited state (thin solid line), as a function of time, is given 
by 

P(t) = \(^(z,t = 0M(z,t))\ 2 . (4) 

This quantity measures excited population in any combination of the first excited states, 
and not only the population in the desired 7r-state. The subsequent Rabi oscillations 
are shown in Fig. 6(a). The corresponding atom number population difference between 
the two wells features 7r-oscillations (Fig. 6(b)) [Hill EH), in which there is almost 
complete exchange of atoms between the two wells (i.e. essentially population exchange 
between the two macroscopically quantum self-trapped states of Fig. 1(c) (iii)). This 
phenomenon acts independently of the linear magnetic field gradient and tends to 
obscure the effect of Josephson counterflow. To demonstrate counterflow in its purest 
form, we thus consider the case (Fig. 6(c)) when the amplitude of such 7r-oscillations 
is suppressed (Fig. 6(d)), such that the majority of the population remains in the left, 
lower well at all times. Note that, in all cases studied here, the coupling with higher 
lying states is found to be negligible. 

The importance of counterflow can be stressed, by first showing typical fractional 
population difference dynamics for a double well condensate in its ground state |27j . 
Application of the external potential at t = creates a potential difference between the 
two wells, with flow occuring towards the lower well (located on the left for A > 0). For 
the relatively slow rates of increase of the potential gradient studied here, the system 
follows the ground eigenstate almost adiabatically (Fig. 7(b)). In contrast, Figs. 7(c)- 
(d) show typical counterflow dynamics induced by the potential gradient A ^ 0. This 
is based on the initial state preparation of Fig. 6(c)-(d) for t < r and subsequent 
free evolution for tq < t < t\ (as shown in Fig. 7(a)), with the potential gradient A 
applied at t = r±. The evolution of the fractional population difference during this 
entire process is shown respectively by the solid lines in Fig. 7(c)- (d) for (c) T\ = lOcuJ 1 
and (d) T\ = 85LUJ 1 . The corresponding times when the potential gradient is applied 
are indicated by open circles in Fig. 7(a),(c)-(d). The effect of Josephson counterflow 
manifests itself clearly in that the population starts being trasferred to the right, higher 
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Figure 6. (a) Initial excited state preparation cycle and subsequent free evolution 
featuring Rabi oscillations (r = 6.5WJ 1 , a — Q.5Jiw±). (b) Corresponding fractional 
population differences, exhibiting 7r-oscillations. (c)-(d) Corresponding plots for 
different initial state preparation cycle (to = Sw^ 1 , a = hu>±). In (d), the 7r-oscillations 
have been suppressed, with macroscopic quantum self-trapping occuring in the left, 
lower well. Other parameters as in Fig. 3(a). 



potential energy well. Note that this is a direct consequence of the imposed magnetic 
field gradient (black curve), and that counterflow here occurs in a direction opposite to, 
and for larger amplitude than, the suppressed 7r-oscillations (grey lines). This indicates 
clearly the distinction between counterflow and 7r-oscillations. 

We now look into typical experimental parameters which allow for demonstration 
of counterflow. In particular, we should investigate whether this effect is observable for 
an experimentally realistic number of atoms in the double-well condensate. The number 
of atoms is given by 

* = l- = T~\f^ ■ W 
47r a Aira y muj± 

It follows that, for given dimensionless nonlinearity g, large condensate atom number can 
be obtained for light, weakly-interacting, transversally weakly-confined systems. Note 
also that the total atom number is independent of the trap aspect ratio, as this cancels 
out for fixed g. Hence, for this effect to be observed clearly, with a large number of atoms, 
one should preferably choose atoms with a small value of a^frn. This will hence yield 
large atom numbers for 7 Li and 23 Na, with the corresponding numbers for 87 Rb, 85 Rb 
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Figure 7. (a) Optimized 7r-state preparation cycle and subsequent free evolution 
featuring Rabi oscillations, (i) Maximum and (ii) minimum achievable efficiency of 
population transfer to first excited state (tq = SuiJ 1 , a = 7ki>j_). (b)-(d) Evolution of 
fractional population difference N as a function of potential gradient A for different 
initial states: (b) Ground state condensate exhibiting usual Josephson flow towards 
lower (left) well, (c) Josephson counterflow (black) after efficient 7r-state preparation, 
induced by application of a potential gradient A = R(t — t\) for t\ = 1 OuiJ 1 (point (i) 
in Fig. (a)). Grey: corresponding evolution in the absence of the potential gradient, 
exhibiting (suppressed) 7r-oscillations. (d) Same as (c), but with t\ — (point 
(ii) in Fig. (a)). Other parameters as in Fig. 3(a). 



considerably smaller. Nonetheless, this effect should be observable for all above species 
in experiments with resolution ability of detecting more than 1000 atoms. For example, 
taking g = An as in Fig. 3(b) and uj± = 2n x 5 Hz, we find the following atom numbers 
in the double well: M = 3300 ( 23 Na) and 1000 ( 87 Rb). This number could be enhanced 
by a factor of 5 when using the nonlinearity of Fig. 3(c) and uj± ~ 2n x 1 Hz, whereas 
further enhancement by a factor of 10 is possible by tuning around a Feshbach resonance 
(e.g. 23 Na, 85 Rb, 133 Cs [45J). In the case of 7 Li, this effect should be observable in a very 
clean manner, since in this case the constraint is placed on the maximum number of 
atoms which can be condensed such that it does not exceed the critical number leading 
to collapse [46J. In the case of 85 Rb, the number of atoms needed to observe Josephson 
counterflow is not likely to exceed the critical value for collapse. The effect of Josephson 
tunnelling on collapse will be investigated in subsequent work |47| . 
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The possibility to demonstrate this effect experimentally will also depend on the 
timescales required for its observation. We consider the case of an applied magnetic 
field gradient R = (10~ 3 — 10~ 2 )(hu±/a±) and a typical transverse trapping frequency 
u± = 2tt x (5 — 100)Hz. This translates into a magnetic field gradient inducing 
a Zeeman shift in the range (1MHz— lGHz)/cm, and a dimensionless timescale of 
cuj 1 = (32 — 1.6)ms. Hence, for the illustrative parameters chosen here, efficient 
preparation of the tt state requires a time To ~ (300 — 150)ms. Observation of Josephson 
counterflow requires monitoring the population difference for at least a further time 
4x P ~ (1-5 s— 75 ms). One notices two competing effects here: For fixed, reasonably 
small, nonlinearity (g < 107r), such that the effect can be clearly observed, one needs 
weak transverse confinement uj± in order to obtain a reasonable number of atoms which 
can be imaged easily. Contrary to this, small u± imply long timescales, such that 
the observation of this effect becomes limited by other factors (e.g. thermal damping 
fHlEHIIEl, 3-body recombination, etc.). The best conditions will hence depend on the 
details of each set-up, with reasonable parameters in the range u± ~ 2-k x (5 — 100)Hz 
and uj z ~ 27r x (1 — 500)Hz, leading to total experimental timescales of t exp ~ 100 ms— 2 
s, and requiring an optimum resolution of a few hundred to a few thousand atoms. 

5. Conclusions 

We have studied the Josephson dynamics of a condensate in a double- well potential in the 
presense of a magnetic field gradient, for a suitably prepared initial antisymmetric state 
featuring equal populations in each well and a 7r-phase slip across the weak link. Under 
appropriate conditions, the atomic current was shown to flow in the direction opposite 
to that of minimum potential energy. This is the opposite behaviour to the 'normal' 
Josephson flow occuring for a system in its ground state. This phenomenon, termed 
here Josephson counterflow, bears close analogies to (metastable) 7r-states observed 
in superfluid- 3 He and (controllable) 7r-junctions in superconducting weak links. We 
have discussed a range of typical parameters for which this effect could be observed in 
atomic Bose-Einstein condensates. Observation of this effect in a 'clean' manner requires 
reasonably light, weakly-interacting atoms under rather weak transverse confinement, 
and techniques to populate the first excited state in a manner such that the 7r-oscillations 
are heavily suppressed. Optimum choice of parameters is an interplay between good 
experimental resolution for detecting few hundred to few thousand atoms, and the 
maximum observation time for which this effect is not affected by other dephasing 
processes. It is important to stress that appearance of this effect is not dependent on 
sign, and only weakly- dependent on strength of the scattering length, applying equally 
well to both attractive and repulsive Bose-Einstein condensates. An alternative suitable 
candidate for observing counterflow might also be found in the recently realized atom 
chips jlHj, which offer excellent control and experimental resolution. Investigation in 
these low dimensional systems should enable observation of Jospehson counterflow, since 
we have found this effect to be only weakly dependent on aspect ratio of the harmonic 
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trap. We believe that the experimental observation of this phenomenon will further 
strengthen the analogies between atomic BEC's and other states exhibiting macroscopic 
phase coherence and controllable Josephson effects. 
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